home *** CD-ROM | disk | FTP | other *** search
/ Enter 2004 January / enter-2004-01.iso / files / maxima-5.9.0.exe / {app} / share / maxima / 5.9.0 / src / numerical / slatec / dasyjy.lisp < prev    next >
Encoding:
Text File  |  2003-02-09  |  51.5 KB  |  1,019 lines
  1. ;;; Compiled by f2cl version 2.0 beta 2002-05-06
  2. ;;; 
  3. ;;; Options: ((:prune-labels nil) (:auto-save t) (:relaxed-array-decls t)
  4. ;;;           (:coerce-assigns :as-needed) (:array-type ':array)
  5. ;;;           (:array-slicing nil) (:declare-common nil)
  6. ;;;           (:float-format double-float))
  7.  
  8. (in-package "SLATEC")
  9.  
  10.  
  11. (let ((tols -6.90775527898214)
  12.       (con1 0.666666666666667)
  13.       (con2 0.33333333333333304)
  14.       (con548 0.10416666666666702)
  15.       (ar (make-array 8 :element-type 'double-float))
  16.       (br (make-array 10 :element-type 'double-float))
  17.       (c (make-array 65 :element-type 'double-float))
  18.       (alfa (make-array 104 :element-type 'double-float))
  19.       (beta (make-array 130 :element-type 'double-float))
  20.       (gama (make-array 26 :element-type 'double-float)))
  21.   (declare (type (array double-float (26)) gama)
  22.            (type (array double-float (130)) beta)
  23.            (type (array double-float (104)) alfa)
  24.            (type (array double-float (65)) c)
  25.            (type (array double-float (10)) br)
  26.            (type (array double-float (8)) ar)
  27.            (type double-float con548 con2 con1 tols))
  28.   (f2cl-lib:fset (f2cl-lib:fref ar (1) ((1 8))) 0.0835503472222222)
  29.   (f2cl-lib:fset (f2cl-lib:fref ar (2) ((1 8))) 0.128226574556327)
  30.   (f2cl-lib:fset (f2cl-lib:fref ar (3) ((1 8))) 0.29184902646414)
  31.   (f2cl-lib:fset (f2cl-lib:fref ar (4) ((1 8))) 0.881627267443758)
  32.   (f2cl-lib:fset (f2cl-lib:fref ar (5) ((1 8))) 3.32140828186277)
  33.   (f2cl-lib:fset (f2cl-lib:fref ar (6) ((1 8))) 14.9957629868626)
  34.   (f2cl-lib:fset (f2cl-lib:fref ar (7) ((1 8))) 78.9230130115865)
  35.   (f2cl-lib:fset (f2cl-lib:fref ar (8) ((1 8))) 474.451538868264)
  36.   (f2cl-lib:fset (f2cl-lib:fref br (1) ((1 10))) -0.145833333333333)
  37.   (f2cl-lib:fset (f2cl-lib:fref br (2) ((1 10))) -0.09874131944444439)
  38.   (f2cl-lib:fset (f2cl-lib:fref br (3) ((1 10))) -0.14331205391589502)
  39.   (f2cl-lib:fset (f2cl-lib:fref br (4) ((1 10))) -0.31722720267841403)
  40.   (f2cl-lib:fset (f2cl-lib:fref br (5) ((1 10))) -0.9424291479571201)
  41.   (f2cl-lib:fset (f2cl-lib:fref br (6) ((1 10))) -3.51120304082635)
  42.   (f2cl-lib:fset (f2cl-lib:fref br (7) ((1 10))) -15.727263620368)
  43.   (f2cl-lib:fset (f2cl-lib:fref br (8) ((1 10))) -82.2814390971859)
  44.   (f2cl-lib:fset (f2cl-lib:fref br (9) ((1 10))) -492.355370523671)
  45.   (f2cl-lib:fset (f2cl-lib:fref br (10) ((1 10))) -3316.2185685479703)
  46.   (f2cl-lib:fset (f2cl-lib:fref c (1) ((1 65))) -0.208333333333333)
  47.   (f2cl-lib:fset (f2cl-lib:fref c (2) ((1 65))) 0.125)
  48.   (f2cl-lib:fset (f2cl-lib:fref c (3) ((1 65))) 0.33420138888888906)
  49.   (f2cl-lib:fset (f2cl-lib:fref c (4) ((1 65))) -0.40104166666666696)
  50.   (f2cl-lib:fset (f2cl-lib:fref c (5) ((1 65))) 0.0703125)
  51.   (f2cl-lib:fset (f2cl-lib:fref c (6) ((1 65))) -1.02581259645062)
  52.   (f2cl-lib:fset (f2cl-lib:fref c (7) ((1 65))) 1.84646267361111)
  53.   (f2cl-lib:fset (f2cl-lib:fref c (8) ((1 65))) -0.8912109375)
  54.   (f2cl-lib:fset (f2cl-lib:fref c (9) ((1 65))) 0.0732421875)
  55.   (f2cl-lib:fset (f2cl-lib:fref c (10) ((1 65))) 4.66958442342625)
  56.   (f2cl-lib:fset (f2cl-lib:fref c (11) ((1 65))) -11.207002616223)
  57.   (f2cl-lib:fset (f2cl-lib:fref c (12) ((1 65))) 8.78912353515625)
  58.   (f2cl-lib:fset (f2cl-lib:fref c (13) ((1 65))) -2.3640869140625)
  59.   (f2cl-lib:fset (f2cl-lib:fref c (14) ((1 65))) 0.112152099609375)
  60.   (f2cl-lib:fset (f2cl-lib:fref c (15) ((1 65))) -28.2120725582002)
  61.   (f2cl-lib:fset (f2cl-lib:fref c (16) ((1 65))) 84.6362176746007)
  62.   (f2cl-lib:fset (f2cl-lib:fref c (17) ((1 65))) -91.81824154323999)
  63.   (f2cl-lib:fset (f2cl-lib:fref c (18) ((1 65))) 42.5349987453885)
  64.   (f2cl-lib:fset (f2cl-lib:fref c (19) ((1 65))) -7.36879435947963)
  65.   (f2cl-lib:fset (f2cl-lib:fref c (20) ((1 65))) 0.22710800170898404)
  66.   (f2cl-lib:fset (f2cl-lib:fref c (21) ((1 65))) 212.57013003921702)
  67.   (f2cl-lib:fset (f2cl-lib:fref c (22) ((1 65))) -765.252468141182)
  68.   (f2cl-lib:fset (f2cl-lib:fref c (23) ((1 65))) 1059.9904525279999)
  69.   (f2cl-lib:fset (f2cl-lib:fref c (24) ((1 65))) -699.5796273761331)
  70.   (f2cl-lib:fset (f2cl-lib:fref c (25) ((1 65))) 218.190511744212)
  71.   (f2cl-lib:fset (f2cl-lib:fref c (26) ((1 65))) -26.4914304869516)
  72.   (f2cl-lib:fset (f2cl-lib:fref c (27) ((1 65))) 0.572501420974731)
  73.   (f2cl-lib:fset (f2cl-lib:fref c (28) ((1 65))) -1919.45766231841)
  74.   (f2cl-lib:fset (f2cl-lib:fref c (29) ((1 65))) 8061.722181737309)
  75.   (f2cl-lib:fset (f2cl-lib:fref c (30) ((1 65))) -13586.5500064341)
  76.   (f2cl-lib:fset (f2cl-lib:fref c (31) ((1 65))) 11655.3933368645)
  77.   (f2cl-lib:fset (f2cl-lib:fref c (32) ((1 65))) -5305.6469786134)
  78.   (f2cl-lib:fset (f2cl-lib:fref c (33) ((1 65))) 1200.90291321635)
  79.   (f2cl-lib:fset (f2cl-lib:fref c (34) ((1 65))) -108.090919788395)
  80.   (f2cl-lib:fset (f2cl-lib:fref c (35) ((1 65))) 1.72772750258446)
  81.   (f2cl-lib:fset (f2cl-lib:fref c (36) ((1 65))) 20204.2913309661)
  82.   (f2cl-lib:fset (f2cl-lib:fref c (37) ((1 65))) -96980.5983886375)
  83.   (f2cl-lib:fset (f2cl-lib:fref c (38) ((1 65))) 192547.001232532)
  84.   (f2cl-lib:fset (f2cl-lib:fref c (39) ((1 65))) -203400.177280416)
  85.   (f2cl-lib:fset (f2cl-lib:fref c (40) ((1 65))) 122200.464983017)
  86.   (f2cl-lib:fset (f2cl-lib:fref c (41) ((1 65))) -41192.6549688976)
  87.   (f2cl-lib:fset (f2cl-lib:fref c (42) ((1 65))) 7109.51430248936)
  88.   (f2cl-lib:fset (f2cl-lib:fref c (43) ((1 65))) -493.915304773088)
  89.   (f2cl-lib:fset (f2cl-lib:fref c (44) ((1 65))) 6.07404200127348)
  90.   (f2cl-lib:fset (f2cl-lib:fref c (45) ((1 65))) -242919.187900551)
  91.   (f2cl-lib:fset (f2cl-lib:fref c (46) ((1 65))) 1311763.6146629802)
  92.   (f2cl-lib:fset (f2cl-lib:fref c (47) ((1 65))) -2998015.9185381103)
  93.   (f2cl-lib:fset (f2cl-lib:fref c (48) ((1 65))) 3763271.2976564)
  94.   (f2cl-lib:fset (f2cl-lib:fref c (49) ((1 65))) -2813563.22658653)
  95.   (f2cl-lib:fset (f2cl-lib:fref c (50) ((1 65))) 1268365.27332162)
  96.   (f2cl-lib:fset (f2cl-lib:fref c (51) ((1 65))) -331645.172484564)
  97.   (f2cl-lib:fset (f2cl-lib:fref c (52) ((1 65))) 45218.7689813627)
  98.   (f2cl-lib:fset (f2cl-lib:fref c (53) ((1 65))) -2499.8304818112097)
  99.   (f2cl-lib:fset (f2cl-lib:fref c (54) ((1 65))) 24.3805296995561)
  100.   (f2cl-lib:fset (f2cl-lib:fref c (55) ((1 65))) 3284469.8530720402)
  101.   (f2cl-lib:fset (f2cl-lib:fref c (56) ((1 65))) -1.9706819118432198e+7)
  102.   (f2cl-lib:fset (f2cl-lib:fref c (57) ((1 65))) 5.09526024926646e+7)
  103.   (f2cl-lib:fset (f2cl-lib:fref c (58) ((1 65))) -7.41051482115327e+7)
  104.   (f2cl-lib:fset (f2cl-lib:fref c (59) ((1 65))) 6.634451227472901e+7)
  105.   (f2cl-lib:fset (f2cl-lib:fref c (60) ((1 65))) -3.7567176660763396e+7)
  106.   (f2cl-lib:fset (f2cl-lib:fref c (61) ((1 65))) 1.32887671664218e+7)
  107.   (f2cl-lib:fset (f2cl-lib:fref c (62) ((1 65))) -2785618.12808645)
  108.   (f2cl-lib:fset (f2cl-lib:fref c (63) ((1 65))) 308186.404612662)
  109.   (f2cl-lib:fset (f2cl-lib:fref c (64) ((1 65))) -13886.089753716999)
  110.   (f2cl-lib:fset (f2cl-lib:fref c (65) ((1 65))) 110.01714026924701)
  111.   (f2cl-lib:fset (f2cl-lib:fref alfa (1 1) ((1 26) (1 4)))
  112.                  -0.00444444444444444)
  113.   (f2cl-lib:fset (f2cl-lib:fref alfa (2 1) ((1 26) (1 4)))
  114.                  -9.220779220779219e-4)
  115.   (f2cl-lib:fset (f2cl-lib:fref alfa (3 1) ((1 26) (1 4)))
  116.                  -8.848928848928851e-5)
  117.   (f2cl-lib:fset (f2cl-lib:fref alfa (4 1) ((1 26) (1 4))) 1.6592768783245e-4)
  118.   (f2cl-lib:fset (f2cl-lib:fref alfa (5 1) ((1 26) (1 4))) 2.46691372741793e-4)
  119.   (f2cl-lib:fset (f2cl-lib:fref alfa (6 1) ((1 26) (1 4))) 2.65995589346255e-4)
  120.   (f2cl-lib:fset (f2cl-lib:fref alfa (7 1) ((1 26) (1 4))) 2.61824297061501e-4)
  121.   (f2cl-lib:fset (f2cl-lib:fref alfa (8 1) ((1 26) (1 4)))
  122.                  2.4873043734465602e-4)
  123.   (f2cl-lib:fset (f2cl-lib:fref alfa (9 1) ((1 26) (1 4)))
  124.                  2.3272104008323197e-4)
  125.   (f2cl-lib:fset (f2cl-lib:fref alfa (10 1) ((1 26) (1 4)))
  126.                  2.1636248571236502e-4)
  127.   (f2cl-lib:fset (f2cl-lib:fref alfa (11 1) ((1 26) (1 4)))
  128.                  2.0073885876275202e-4)
  129.   (f2cl-lib:fset (f2cl-lib:fref alfa (12 1) ((1 26) (1 4)))
  130.                  1.8626763663754503e-4)
  131.   (f2cl-lib:fset (f2cl-lib:fref alfa (13 1) ((1 26) (1 4)))
  132.                  1.73060775917876e-4)
  133.   (f2cl-lib:fset (f2cl-lib:fref alfa (14 1) ((1 26) (1 4)))
  134.                  1.61091705929016e-4)
  135.   (f2cl-lib:fset (f2cl-lib:fref alfa (15 1) ((1 26) (1 4)))
  136.                  1.50274774160908e-4)
  137.   (f2cl-lib:fset (f2cl-lib:fref alfa (16 1) ((1 26) (1 4)))
  138.                  1.4050349739127002e-4)
  139.   (f2cl-lib:fset (f2cl-lib:fref alfa (17 1) ((1 26) (1 4)))
  140.                  1.3166881654592302e-4)
  141.   (f2cl-lib:fset (f2cl-lib:fref alfa (18 1) ((1 26) (1 4)))
  142.                  1.23667445598253e-4)
  143.   (f2cl-lib:fset (f2cl-lib:fref alfa (19 1) ((1 26) (1 4)))
  144.                  1.1640527147473803e-4)
  145.   (f2cl-lib:fset (f2cl-lib:fref alfa (20 1) ((1 26) (1 4)))
  146.                  1.09798298372713e-4)
  147.   (f2cl-lib:fset (f2cl-lib:fref alfa (21 1) ((1 26) (1 4)))
  148.                  1.0377241042299298e-4)
  149.   (f2cl-lib:fset (f2cl-lib:fref alfa (22 1) ((1 26) (1 4)))
  150.                  9.826260783693632e-5)
  151.   (f2cl-lib:fset (f2cl-lib:fref alfa (23 1) ((1 26) (1 4)))
  152.                  9.321205172495033e-5)
  153.   (f2cl-lib:fset (f2cl-lib:fref alfa (24 1) ((1 26) (1 4)))
  154.                  8.857108524787122e-5)
  155.   (f2cl-lib:fset (f2cl-lib:fref alfa (25 1) ((1 26) (1 4))) 8.429631057157e-5)
  156.   (f2cl-lib:fset (f2cl-lib:fref alfa (26 1) ((1 26) (1 4)))
  157.                  8.03497548407791e-5)
  158.   (f2cl-lib:fset (f2cl-lib:fref alfa (1 2) ((1 26) (1 4)))
  159.                  6.937355413545891e-4)
  160.   (f2cl-lib:fset (f2cl-lib:fref alfa (2 2) ((1 26) (1 4))) 2.32241745182922e-4)
  161.   (f2cl-lib:fset (f2cl-lib:fref alfa (3 2) ((1 26) (1 4)))
  162.                  -1.41986273556691e-5)
  163.   (f2cl-lib:fset (f2cl-lib:fref alfa (4 2) ((1 26) (1 4)))
  164.                  -1.1644493167204902e-4)
  165.   (f2cl-lib:fset (f2cl-lib:fref alfa (5 2) ((1 26) (1 4)))
  166.                  -1.5080355805304899e-4)
  167.   (f2cl-lib:fset (f2cl-lib:fref alfa (6 2) ((1 26) (1 4)))
  168.                  -1.55121924918096e-4)
  169.   (f2cl-lib:fset (f2cl-lib:fref alfa (7 2) ((1 26) (1 4)))
  170.                  -1.46809756646466e-4)
  171.   (f2cl-lib:fset (f2cl-lib:fref alfa (8 2) ((1 26) (1 4)))
  172.                  -1.3381550386749097e-4)
  173.   (f2cl-lib:fset (f2cl-lib:fref alfa (9 2) ((1 26) (1 4)))
  174.                  -1.1974497568425399e-4)
  175.   (f2cl-lib:fset (f2cl-lib:fref alfa (10 2) ((1 26) (1 4)))
  176.                  -1.06184319207974e-4)
  177.   (f2cl-lib:fset (f2cl-lib:fref alfa (11 2) ((1 26) (1 4)))
  178.                  -9.376995498911941e-5)
  179.   (f2cl-lib:fset (f2cl-lib:fref alfa (12 2) ((1 26) (1 4)))
  180.                  -8.26923045588193e-5)
  181.   (f2cl-lib:fset (f2cl-lib:fref alfa (13 2) ((1 26) (1 4)))
  182.                  -7.29374348155221e-5)
  183.   (f2cl-lib:fset (f2cl-lib:fref alfa (14 2) ((1 26) (1 4)))
  184.                  -6.44042357721016e-5)
  185.   (f2cl-lib:fset (f2cl-lib:fref alfa (15 2) ((1 26) (1 4)))
  186.                  -5.69611566009369e-5)
  187.   (f2cl-lib:fset (f2cl-lib:fref alfa (16 2) ((1 26) (1 4)))
  188.                  -5.047310443035621e-5)
  189.   (f2cl-lib:fset (f2cl-lib:fref alfa (17 2) ((1 26) (1 4)))
  190.                  -4.481348680088831e-5)
  191.   (f2cl-lib:fset (f2cl-lib:fref alfa (18 2) ((1 26) (1 4)))
  192.                  -3.9868872771759906e-5)
  193.   (f2cl-lib:fset (f2cl-lib:fref alfa (19 2) ((1 26) (1 4)))
  194.                  -3.55400532972042e-5)
  195.   (f2cl-lib:fset (f2cl-lib:fref alfa (20 2) ((1 26) (1 4)))
  196.                  -3.17414256609022e-5)
  197.   (f2cl-lib:fset (f2cl-lib:fref alfa (21 2) ((1 26) (1 4)))
  198.                  -2.83996793904175e-5)
  199.   (f2cl-lib:fset (f2cl-lib:fref alfa (22 2) ((1 26) (1 4)))
  200.                  -2.54522720634871e-5)
  201.   (f2cl-lib:fset (f2cl-lib:fref alfa (23 2) ((1 26) (1 4)))
  202.                  -2.2845929716472502e-5)
  203.   (f2cl-lib:fset (f2cl-lib:fref alfa (24 2) ((1 26) (1 4)))
  204.                  -2.05352753106481e-5)
  205.   (f2cl-lib:fset (f2cl-lib:fref alfa (25 2) ((1 26) (1 4)))
  206.                  -1.8481621762766604e-5)
  207.   (f2cl-lib:fset (f2cl-lib:fref alfa (26 2) ((1 26) (1 4)))
  208.                  -1.6651933002139402e-5)
  209.   (f2cl-lib:fset (f2cl-lib:fref alfa (1 3) ((1 26) (1 4)))
  210.                  -3.54211971457744e-4)
  211.   (f2cl-lib:fset (f2cl-lib:fref alfa (2 3) ((1 26) (1 4)))
  212.                  -1.56161263945159e-4)
  213.   (f2cl-lib:fset (f2cl-lib:fref alfa (3 3) ((1 26) (1 4)))
  214.                  3.0446550359493607e-5)
  215.   (f2cl-lib:fset (f2cl-lib:fref alfa (4 3) ((1 26) (1 4)))
  216.                  1.3019865577324302e-4)
  217.   (f2cl-lib:fset (f2cl-lib:fref alfa (5 3) ((1 26) (1 4)))
  218.                  1.6747110669971202e-4)
  219.   (f2cl-lib:fset (f2cl-lib:fref alfa (6 3) ((1 26) (1 4)))
  220.                  1.7022258768359302e-4)
  221.   (f2cl-lib:fset (f2cl-lib:fref alfa (7 3) ((1 26) (1 4))) 1.56501427608595e-4)
  222.   (f2cl-lib:fset (f2cl-lib:fref alfa (8 3) ((1 26) (1 4))) 1.36339170977445e-4)
  223.   (f2cl-lib:fset (f2cl-lib:fref alfa (9 3) ((1 26) (1 4))) 1.14886692029825e-4)
  224.   (f2cl-lib:fset (f2cl-lib:fref alfa (10 3) ((1 26) (1 4)))
  225.                  9.45869093034688e-5)
  226.   (f2cl-lib:fset (f2cl-lib:fref alfa (11 3) ((1 26) (1 4)))
  227.                  7.64498419250898e-5)
  228.   (f2cl-lib:fset (f2cl-lib:fref alfa (12 3) ((1 26) (1 4)))
  229.                  6.075703349651971e-5)
  230.   (f2cl-lib:fset (f2cl-lib:fref alfa (13 3) ((1 26) (1 4)))
  231.                  4.743942992905091e-5)
  232.   (f2cl-lib:fset (f2cl-lib:fref alfa (14 3) ((1 26) (1 4)))
  233.                  3.62757512005344e-5)
  234.   (f2cl-lib:fset (f2cl-lib:fref alfa (15 3) ((1 26) (1 4)))
  235.                  2.69939714979225e-5)
  236.   (f2cl-lib:fset (f2cl-lib:fref alfa (16 3) ((1 26) (1 4)))
  237.                  1.93210938247939e-5)
  238.   (f2cl-lib:fset (f2cl-lib:fref alfa (17 3) ((1 26) (1 4)))
  239.                  1.3005667479396302e-5)
  240.   (f2cl-lib:fset (f2cl-lib:fref alfa (18 3) ((1 26) (1 4)))
  241.                  7.826208667444968e-6)
  242.   (f2cl-lib:fset (f2cl-lib:fref alfa (19 3) ((1 26) (1 4)))
  243.                  3.59257485819352e-6)
  244.   (f2cl-lib:fset (f2cl-lib:fref alfa (20 3) ((1 26) (1 4)))
  245.                  1.4404004981425203e-7)
  246.   (f2cl-lib:fset (f2cl-lib:fref alfa (21 3) ((1 26) (1 4)))
  247.                  -2.65396769697939e-6)
  248.   (f2cl-lib:fset (f2cl-lib:fref alfa (22 3) ((1 26) (1 4)))
  249.                  -4.913468670984859e-6)
  250.   (f2cl-lib:fset (f2cl-lib:fref alfa (23 3) ((1 26) (1 4)))
  251.                  -6.72739296091248e-6)
  252.   (f2cl-lib:fset (f2cl-lib:fref alfa (24 3) ((1 26) (1 4)))
  253.                  -8.17269379678658e-6)
  254.   (f2cl-lib:fset (f2cl-lib:fref alfa (25 3) ((1 26) (1 4)))
  255.                  -9.313047150935608e-6)
  256.   (f2cl-lib:fset (f2cl-lib:fref alfa (26 3) ((1 26) (1 4)))
  257.                  -1.0201141879801601e-5)
  258.   (f2cl-lib:fset (f2cl-lib:fref alfa (1 4) ((1 26) (1 4))) 3.78194199201773e-4)
  259.   (f2cl-lib:fset (f2cl-lib:fref alfa (2 4) ((1 26) (1 4))) 2.02471952761816e-4)
  260.   (f2cl-lib:fset (f2cl-lib:fref alfa (3 4) ((1 26) (1 4)))
  261.                  -6.379385063188621e-5)
  262.   (f2cl-lib:fset (f2cl-lib:fref alfa (4 4) ((1 26) (1 4)))
  263.                  -2.3859823060300597e-4)
  264.   (f2cl-lib:fset (f2cl-lib:fref alfa (5 4) ((1 26) (1 4)))
  265.                  -3.1091625602736206e-4)
  266.   (f2cl-lib:fset (f2cl-lib:fref alfa (6 4) ((1 26) (1 4)))
  267.                  -3.13680115247576e-4)
  268.   (f2cl-lib:fset (f2cl-lib:fref alfa (7 4) ((1 26) (1 4)))
  269.                  -2.78950273791323e-4)
  270.   (f2cl-lib:fset (f2cl-lib:fref alfa (8 4) ((1 26) (1 4)))
  271.                  -2.28564082619141e-4)
  272.   (f2cl-lib:fset (f2cl-lib:fref alfa (9 4) ((1 26) (1 4)))
  273.                  -1.7524528034084702e-4)
  274.   (f2cl-lib:fset (f2cl-lib:fref alfa (10 4) ((1 26) (1 4)))
  275.                  -1.2554406306069002e-4)
  276.   (f2cl-lib:fset (f2cl-lib:fref alfa (11 4) ((1 26) (1 4)))
  277.                  -8.229828728202081e-5)
  278.   (f2cl-lib:fset (f2cl-lib:fref alfa (12 4) ((1 26) (1 4)))
  279.                  -4.62860730588116e-5)
  280.   (f2cl-lib:fset (f2cl-lib:fref alfa (13 4) ((1 26) (1 4)))
  281.                  -1.7233430236696204e-5)
  282.   (f2cl-lib:fset (f2cl-lib:fref alfa (14 4) ((1 26) (1 4)))
  283.                  5.606904823046021e-6)
  284.   (f2cl-lib:fset (f2cl-lib:fref alfa (15 4) ((1 26) (1 4)))
  285.                  2.3139544314828706e-5)
  286.   (f2cl-lib:fset (f2cl-lib:fref alfa (16 4) ((1 26) (1 4)))
  287.                  3.6264274585679407e-5)
  288.   (f2cl-lib:fset (f2cl-lib:fref alfa (17 4) ((1 26) (1 4)))
  289.                  4.5800612449018896e-5)
  290.   (f2cl-lib:fset (f2cl-lib:fref alfa (18 4) ((1 26) (1 4)))
  291.                  5.245952949591141e-5)
  292.   (f2cl-lib:fset (f2cl-lib:fref alfa (19 4) ((1 26) (1 4)))
  293.                  5.683962085458151e-5)
  294.   (f2cl-lib:fset (f2cl-lib:fref alfa (20 4) ((1 26) (1 4)))
  295.                  5.94349820393104e-5)
  296.   (f2cl-lib:fset (f2cl-lib:fref alfa (21 4) ((1 26) (1 4)))
  297.                  6.06478527578422e-5)
  298.   (f2cl-lib:fset (f2cl-lib:fref alfa (22 4) ((1 26) (1 4)))
  299.                  6.080239077884361e-5)
  300.   (f2cl-lib:fset (f2cl-lib:fref alfa (23 4) ((1 26) (1 4))) 6.0157789453946e-5)
  301.   (f2cl-lib:fset (f2cl-lib:fref alfa (24 4) ((1 26) (1 4)))
  302.                  5.89199657344698e-5)
  303.   (f2cl-lib:fset (f2cl-lib:fref alfa (25 4) ((1 26) (1 4)))
  304.                  5.72515823777593e-5)
  305.   (f2cl-lib:fset (f2cl-lib:fref alfa (26 4) ((1 26) (1 4)))
  306.                  5.52804375585853e-5)
  307.   (f2cl-lib:fset (f2cl-lib:fref beta (1 1) ((1 26) (1 5))) 0.0179988721413553)
  308.   (f2cl-lib:fset (f2cl-lib:fref beta (2 1) ((1 26) (1 5))) 0.00559964911064388)
  309.   (f2cl-lib:fset (f2cl-lib:fref beta (3 1) ((1 26) (1 5)))
  310.                  0.0028850140223113303)
  311.   (f2cl-lib:fset (f2cl-lib:fref beta (4 1) ((1 26) (1 5))) 0.00180096606761054)
  312.   (f2cl-lib:fset (f2cl-lib:fref beta (5 1) ((1 26) (1 5)))
  313.                  0.0012475311058919901)
  314.   (f2cl-lib:fset (f2cl-lib:fref beta (6 1) ((1 26) (1 5)))
  315.                  9.228788765729381e-4)
  316.   (f2cl-lib:fset (f2cl-lib:fref beta (7 1) ((1 26) (1 5))) 7.14430421727287e-4)
  317.   (f2cl-lib:fset (f2cl-lib:fref beta (8 1) ((1 26) (1 5)))
  318.                  5.717872817897051e-4)
  319.   (f2cl-lib:fset (f2cl-lib:fref beta (9 1) ((1 26) (1 5)))
  320.                  4.6943100760648204e-4)
  321.   (f2cl-lib:fset (f2cl-lib:fref beta (10 1) ((1 26) (1 5)))
  322.                  3.9323283546291704e-4)
  323.   (f2cl-lib:fset (f2cl-lib:fref beta (11 1) ((1 26) (1 5)))
  324.                  3.34818889318298e-4)
  325.   (f2cl-lib:fset (f2cl-lib:fref beta (12 1) ((1 26) (1 5)))
  326.                  2.88952148495752e-4)
  327.   (f2cl-lib:fset (f2cl-lib:fref beta (13 1) ((1 26) (1 5)))
  328.                  2.5221161554957305e-4)
  329.   (f2cl-lib:fset (f2cl-lib:fref beta (14 1) ((1 26) (1 5)))
  330.                  2.2228058079888302e-4)
  331.   (f2cl-lib:fset (f2cl-lib:fref beta (15 1) ((1 26) (1 5)))
  332.                  1.97541838033063e-4)
  333.   (f2cl-lib:fset (f2cl-lib:fref beta (16 1) ((1 26) (1 5)))
  334.                  1.7683685501971802e-4)
  335.   (f2cl-lib:fset (f2cl-lib:fref beta (17 1) ((1 26) (1 5)))
  336.                  1.59316899661821e-4)
  337.   (f2cl-lib:fset (f2cl-lib:fref beta (18 1) ((1 26) (1 5)))
  338.                  1.44347930197334e-4)
  339.   (f2cl-lib:fset (f2cl-lib:fref beta (19 1) ((1 26) (1 5)))
  340.                  1.3144806811996498e-4)
  341.   (f2cl-lib:fset (f2cl-lib:fref beta (20 1) ((1 26) (1 5)))
  342.                  1.2024544494930303e-4)
  343.   (f2cl-lib:fset (f2cl-lib:fref beta (21 1) ((1 26) (1 5)))
  344.                  1.10449144504599e-4)
  345.   (f2cl-lib:fset (f2cl-lib:fref beta (22 1) ((1 26) (1 5)))
  346.                  1.0182877074056702e-4)
  347.   (f2cl-lib:fset (f2cl-lib:fref beta (23 1) ((1 26) (1 5)))
  348.                  9.419982242042381e-5)
  349.   (f2cl-lib:fset (f2cl-lib:fref beta (24 1) ((1 26) (1 5)))
  350.                  8.741305457538341e-5)
  351.   (f2cl-lib:fset (f2cl-lib:fref beta (25 1) ((1 26) (1 5)))
  352.                  8.13466262162801e-5)
  353.   (f2cl-lib:fset (f2cl-lib:fref beta (26 1) ((1 26) (1 5)))
  354.                  7.590022696462191e-5)
  355.   (f2cl-lib:fset (f2cl-lib:fref beta (1 2) ((1 26) (1 5)))
  356.                  -0.00149282953213429)
  357.   (f2cl-lib:fset (f2cl-lib:fref beta (2 2) ((1 26) (1 5)))
  358.                  -8.78204709546389e-4)
  359.   (f2cl-lib:fset (f2cl-lib:fref beta (3 2) ((1 26) (1 5)))
  360.                  -5.02916549572035e-4)
  361.   (f2cl-lib:fset (f2cl-lib:fref beta (4 2) ((1 26) (1 5)))
  362.                  -2.94822138512746e-4)
  363.   (f2cl-lib:fset (f2cl-lib:fref beta (5 2) ((1 26) (1 5)))
  364.                  -1.7546399697078302e-4)
  365.   (f2cl-lib:fset (f2cl-lib:fref beta (6 2) ((1 26) (1 5)))
  366.                  -1.0400855046081603e-4)
  367.   (f2cl-lib:fset (f2cl-lib:fref beta (7 2) ((1 26) (1 5)))
  368.                  -5.96141953046458e-5)
  369.   (f2cl-lib:fset (f2cl-lib:fref beta (8 2) ((1 26) (1 5)))
  370.                  -3.12038929076098e-5)
  371.   (f2cl-lib:fset (f2cl-lib:fref beta (9 2) ((1 26) (1 5)))
  372.                  -1.2608973598023004e-5)
  373.   (f2cl-lib:fset (f2cl-lib:fref beta (10 2) ((1 26) (1 5)))
  374.                  -2.4289260857573003e-7)
  375.   (f2cl-lib:fset (f2cl-lib:fref beta (11 2) ((1 26) (1 5)))
  376.                  8.059961654142739e-6)
  377.   (f2cl-lib:fset (f2cl-lib:fref beta (12 2) ((1 26) (1 5)))
  378.                  1.3650700926214704e-5)
  379.   (f2cl-lib:fset (f2cl-lib:fref beta (13 2) ((1 26) (1 5)))
  380.                  1.73964125472926e-5)
  381.   (f2cl-lib:fset (f2cl-lib:fref beta (14 2) ((1 26) (1 5)))
  382.                  1.9867297884213403e-5)
  383.   (f2cl-lib:fset (f2cl-lib:fref beta (15 2) ((1 26) (1 5)))
  384.                  2.1446326379082298e-5)
  385.   (f2cl-lib:fset (f2cl-lib:fref beta (16 2) ((1 26) (1 5)))
  386.                  2.2395465923245705e-5)
  387.   (f2cl-lib:fset (f2cl-lib:fref beta (17 2) ((1 26) (1 5)))
  388.                  2.2896778381471306e-5)
  389.   (f2cl-lib:fset (f2cl-lib:fref beta (18 2) ((1 26) (1 5)))
  390.                  2.3078538981117808e-5)
  391.   (f2cl-lib:fset (f2cl-lib:fref beta (19 2) ((1 26) (1 5)))
  392.                  2.3032197608090899e-5)
  393.   (f2cl-lib:fset (f2cl-lib:fref beta (20 2) ((1 26) (1 5)))
  394.                  2.28236073720349e-5)
  395.   (f2cl-lib:fset (f2cl-lib:fref beta (21 2) ((1 26) (1 5)))
  396.                  2.2500588110529204e-5)
  397.   (f2cl-lib:fset (f2cl-lib:fref beta (22 2) ((1 26) (1 5)))
  398.                  2.2098101536199102e-5)
  399.   (f2cl-lib:fset (f2cl-lib:fref beta (23 2) ((1 26) (1 5)))
  400.                  2.1641842744810397e-5)
  401.   (f2cl-lib:fset (f2cl-lib:fref beta (24 2) ((1 26) (1 5)))
  402.                  2.1150764925622098e-5)
  403.   (f2cl-lib:fset (f2cl-lib:fref beta (25 2) ((1 26) (1 5)))
  404.                  2.0638874978217098e-5)
  405.   (f2cl-lib:fset (f2cl-lib:fref beta (26 2) ((1 26) (1 5)))
  406.                  2.0116524199708202e-5)
  407.   (f2cl-lib:fset (f2cl-lib:fref beta (1 3) ((1 26) (1 5))) 5.52213076721293e-4)
  408.   (f2cl-lib:fset (f2cl-lib:fref beta (2 3) ((1 26) (1 5))) 4.47932581552385e-4)
  409.   (f2cl-lib:fset (f2cl-lib:fref beta (3 3) ((1 26) (1 5)))
  410.                  2.7952065399202103e-4)
  411.   (f2cl-lib:fset (f2cl-lib:fref beta (4 3) ((1 26) (1 5))) 1.52468156198447e-4)
  412.   (f2cl-lib:fset (f2cl-lib:fref beta (5 3) ((1 26) (1 5))) 6.93271105657044e-5)
  413.   (f2cl-lib:fset (f2cl-lib:fref beta (6 3) ((1 26) (1 5))) 1.76258683069991e-5)
  414.   (f2cl-lib:fset (f2cl-lib:fref beta (7 3) ((1 26) (1 5)))
  415.                  -1.35744996343269e-5)
  416.   (f2cl-lib:fset (f2cl-lib:fref beta (8 3) ((1 26) (1 5)))
  417.                  -3.1797241335042703e-5)
  418.   (f2cl-lib:fset (f2cl-lib:fref beta (9 3) ((1 26) (1 5)))
  419.                  -4.1886186169669304e-5)
  420.   (f2cl-lib:fset (f2cl-lib:fref beta (10 3) ((1 26) (1 5)))
  421.                  -4.6900488937914103e-5)
  422.   (f2cl-lib:fset (f2cl-lib:fref beta (11 3) ((1 26) (1 5)))
  423.                  -4.87665447413787e-5)
  424.   (f2cl-lib:fset (f2cl-lib:fref beta (12 3) ((1 26) (1 5)))
  425.                  -4.87010031186735e-5)
  426.   (f2cl-lib:fset (f2cl-lib:fref beta (13 3) ((1 26) (1 5)))
  427.                  -4.7475562089008705e-5)
  428.   (f2cl-lib:fset (f2cl-lib:fref beta (14 3) ((1 26) (1 5)))
  429.                  -4.5581305813862805e-5)
  430.   (f2cl-lib:fset (f2cl-lib:fref beta (15 3) ((1 26) (1 5)))
  431.                  -4.33309644511266e-5)
  432.   (f2cl-lib:fset (f2cl-lib:fref beta (16 3) ((1 26) (1 5)))
  433.                  -4.0923019315775006e-5)
  434.   (f2cl-lib:fset (f2cl-lib:fref beta (17 3) ((1 26) (1 5)))
  435.                  -3.84822638603221e-5)
  436.   (f2cl-lib:fset (f2cl-lib:fref beta (18 3) ((1 26) (1 5)))
  437.                  -3.60857167535411e-5)
  438.   (f2cl-lib:fset (f2cl-lib:fref beta (19 3) ((1 26) (1 5)))
  439.                  -3.3779330612336705e-5)
  440.   (f2cl-lib:fset (f2cl-lib:fref beta (20 3) ((1 26) (1 5)))
  441.                  -3.1588856077211e-5)
  442.   (f2cl-lib:fset (f2cl-lib:fref beta (21 3) ((1 26) (1 5)))
  443.                  -2.9526956175080704e-5)
  444.   (f2cl-lib:fset (f2cl-lib:fref beta (22 3) ((1 26) (1 5)))
  445.                  -2.7597891482833603e-5)
  446.   (f2cl-lib:fset (f2cl-lib:fref beta (23 3) ((1 26) (1 5)))
  447.                  -2.58006174666884e-5)
  448.   (f2cl-lib:fset (f2cl-lib:fref beta (24 3) ((1 26) (1 5)))
  449.                  -2.4130835676128e-5)
  450.   (f2cl-lib:fset (f2cl-lib:fref beta (25 3) ((1 26) (1 5)))
  451.                  -2.2582350951834607e-5)
  452.   (f2cl-lib:fset (f2cl-lib:fref beta (26 3) ((1 26) (1 5)))
  453.                  -2.11479656768913e-5)
  454.   (f2cl-lib:fset (f2cl-lib:fref beta (1 4) ((1 26) (1 5))) -4.7461779655996e-4)
  455.   (f2cl-lib:fset (f2cl-lib:fref beta (2 4) ((1 26) (1 5)))
  456.                  -4.778645671473211e-4)
  457.   (f2cl-lib:fset (f2cl-lib:fref beta (3 4) ((1 26) (1 5)))
  458.                  -3.20390228067038e-4)
  459.   (f2cl-lib:fset (f2cl-lib:fref beta (4 4) ((1 26) (1 5)))
  460.                  -1.6110501611996203e-4)
  461.   (f2cl-lib:fset (f2cl-lib:fref beta (5 4) ((1 26) (1 5)))
  462.                  -4.2577810128543503e-5)
  463.   (f2cl-lib:fset (f2cl-lib:fref beta (6 4) ((1 26) (1 5))) 3.44571294294968e-5)
  464.   (f2cl-lib:fset (f2cl-lib:fref beta (7 4) ((1 26) (1 5)))
  465.                  7.970926840756751e-5)
  466.   (f2cl-lib:fset (f2cl-lib:fref beta (8 4) ((1 26) (1 5))) 1.03138236708272e-4)
  467.   (f2cl-lib:fset (f2cl-lib:fref beta (9 4) ((1 26) (1 5)))
  468.                  1.1246677526220401e-4)
  469.   (f2cl-lib:fset (f2cl-lib:fref beta (10 4) ((1 26) (1 5)))
  470.                  1.1310364210848103e-4)
  471.   (f2cl-lib:fset (f2cl-lib:fref beta (11 4) ((1 26) (1 5)))
  472.                  1.0865163484877402e-4)
  473.   (f2cl-lib:fset (f2cl-lib:fref beta (12 4) ((1 26) (1 5)))
  474.                  1.01437951597662e-4)
  475.   (f2cl-lib:fset (f2cl-lib:fref beta (13 4) ((1 26) (1 5)))
  476.                  9.292983965933641e-5)
  477.   (f2cl-lib:fset (f2cl-lib:fref beta (14 4) ((1 26) (1 5)))
  478.                  8.402931330160901e-5)
  479.   (f2cl-lib:fset (f2cl-lib:fref beta (15 4) ((1 26) (1 5)))
  480.                  7.527279913491342e-5)
  481.   (f2cl-lib:fset (f2cl-lib:fref beta (16 4) ((1 26) (1 5)))
  482.                  6.696325219757311e-5)
  483.   (f2cl-lib:fset (f2cl-lib:fref beta (17 4) ((1 26) (1 5)))
  484.                  5.925645473231951e-5)
  485.   (f2cl-lib:fset (f2cl-lib:fref beta (18 4) ((1 26) (1 5)))
  486.                  5.22169308826976e-5)
  487.   (f2cl-lib:fset (f2cl-lib:fref beta (19 4) ((1 26) (1 5)))
  488.                  4.58539485165361e-5)
  489.   (f2cl-lib:fset (f2cl-lib:fref beta (20 4) ((1 26) (1 5)))
  490.                  4.01445513891487e-5)
  491.   (f2cl-lib:fset (f2cl-lib:fref beta (21 4) ((1 26) (1 5)))
  492.                  3.50481730031328e-5)
  493.   (f2cl-lib:fset (f2cl-lib:fref beta (22 4) ((1 26) (1 5)))
  494.                  3.0515799503434704e-5)
  495.   (f2cl-lib:fset (f2cl-lib:fref beta (23 4) ((1 26) (1 5)))
  496.                  2.64956119950516e-5)
  497.   (f2cl-lib:fset (f2cl-lib:fref beta (24 4) ((1 26) (1 5)))
  498.                  2.2936363369099802e-5)
  499.   (f2cl-lib:fset (f2cl-lib:fref beta (25 4) ((1 26) (1 5)))
  500.                  1.97893056664022e-5)
  501.   (f2cl-lib:fset (f2cl-lib:fref beta (26 4) ((1 26) (1 5)))
  502.                  1.70091984636413e-5)
  503.   (f2cl-lib:fset (f2cl-lib:fref beta (1 5) ((1 26) (1 5))) 7.36465810572578e-4)
  504.   (f2cl-lib:fset (f2cl-lib:fref beta (2 5) ((1 26) (1 5))) 8.72790805146194e-4)
  505.   (f2cl-lib:fset (f2cl-lib:fref beta (3 5) ((1 26) (1 5)))
  506.                  6.226148625731349e-4)
  507.   (f2cl-lib:fset (f2cl-lib:fref beta (4 5) ((1 26) (1 5))) 2.85998154194304e-4)
  508.   (f2cl-lib:fset (f2cl-lib:fref beta (5 5) ((1 26) (1 5))) 3.84737672879366e-6)
  509.   (f2cl-lib:fset (f2cl-lib:fref beta (6 5) ((1 26) (1 5)))
  510.                  -1.8790600363697202e-4)
  511.   (f2cl-lib:fset (f2cl-lib:fref beta (7 5) ((1 26) (1 5)))
  512.                  -2.97603646594555e-4)
  513.   (f2cl-lib:fset (f2cl-lib:fref beta (8 5) ((1 26) (1 5)))
  514.                  -3.45998126832656e-4)
  515.   (f2cl-lib:fset (f2cl-lib:fref beta (9 5) ((1 26) (1 5)))
  516.                  -3.53382470916038e-4)
  517.   (f2cl-lib:fset (f2cl-lib:fref beta (10 5) ((1 26) (1 5)))
  518.                  -3.3571563577504904e-4)
  519.   (f2cl-lib:fset (f2cl-lib:fref beta (11 5) ((1 26) (1 5)))
  520.                  -3.0432112478904e-4)
  521.   (f2cl-lib:fset (f2cl-lib:fref beta (12 5) ((1 26) (1 5)))
  522.                  -2.667227230476131e-4)
  523.   (f2cl-lib:fset (f2cl-lib:fref beta (13 5) ((1 26) (1 5)))
  524.                  -2.2765421412282e-4)
  525.   (f2cl-lib:fset (f2cl-lib:fref beta (14 5) ((1 26) (1 5)))
  526.                  -1.8992261185456202e-4)
  527.   (f2cl-lib:fset (f2cl-lib:fref beta (15 5) ((1 26) (1 5)))
  528.                  -1.55058918599094e-4)
  529.   (f2cl-lib:fset (f2cl-lib:fref beta (16 5) ((1 26) (1 5)))
  530.                  -1.2377824076187402e-4)
  531.   (f2cl-lib:fset (f2cl-lib:fref beta (17 5) ((1 26) (1 5)))
  532.                  -9.629261477176442e-5)
  533.   (f2cl-lib:fset (f2cl-lib:fref beta (18 5) ((1 26) (1 5)))
  534.                  -7.251783277144251e-5)
  535.   (f2cl-lib:fset (f2cl-lib:fref beta (19 5) ((1 26) (1 5)))
  536.                  -5.22070028895634e-5)
  537.   (f2cl-lib:fset (f2cl-lib:fref beta (20 5) ((1 26) (1 5)))
  538.                  -3.50347750511901e-5)
  539.   (f2cl-lib:fset (f2cl-lib:fref beta (21 5) ((1 26) (1 5)))
  540.                  -2.06489761035552e-5)
  541.   (f2cl-lib:fset (f2cl-lib:fref beta (22 5) ((1 26) (1 5)))
  542.                  -8.701060968497669e-6)
  543.   (f2cl-lib:fset (f2cl-lib:fref beta (23 5) ((1 26) (1 5)))
  544.                  1.1369868667509997e-6)
  545.   (f2cl-lib:fset (f2cl-lib:fref beta (24 5) ((1 26) (1 5)))
  546.                  9.16426474122779e-6)
  547.   (f2cl-lib:fset (f2cl-lib:fref beta (25 5) ((1 26) (1 5)))
  548.                  1.56477785428873e-5)
  549.   (f2cl-lib:fset (f2cl-lib:fref beta (26 5) ((1 26) (1 5)))
  550.                  2.0822362948246703e-5)
  551.   (f2cl-lib:fset (f2cl-lib:fref gama (1) ((1 26))) 0.629960524947437)
  552.   (f2cl-lib:fset (f2cl-lib:fref gama (2) ((1 26))) 0.251984209978975)
  553.   (f2cl-lib:fset (f2cl-lib:fref gama (3) ((1 26))) 0.154790300415656)
  554.   (f2cl-lib:fset (f2cl-lib:fref gama (4) ((1 26))) 0.110713062416159)
  555.   (f2cl-lib:fset (f2cl-lib:fref gama (5) ((1 26))) 0.08573093955273951)
  556.   (f2cl-lib:fset (f2cl-lib:fref gama (6) ((1 26))) 0.0697161316958684)
  557.   (f2cl-lib:fset (f2cl-lib:fref gama (7) ((1 26))) 0.0586085671893714)
  558.   (f2cl-lib:fset (f2cl-lib:fref gama (8) ((1 26))) 0.0504698873536311)
  559.   (f2cl-lib:fset (f2cl-lib:fref gama (9) ((1 26))) 0.0442600580689155)
  560.   (f2cl-lib:fset (f2cl-lib:fref gama (10) ((1 26))) 0.039372066154351)
  561.   (f2cl-lib:fset (f2cl-lib:fref gama (11) ((1 26))) 0.0354283195924455)
  562.   (f2cl-lib:fset (f2cl-lib:fref gama (12) ((1 26))) 0.032181885750209804)
  563.   (f2cl-lib:fset (f2cl-lib:fref gama (13) ((1 26))) 0.0294646240791158)
  564.   (f2cl-lib:fset (f2cl-lib:fref gama (14) ((1 26))) 0.027158167711293403)
  565.   (f2cl-lib:fset (f2cl-lib:fref gama (15) ((1 26))) 0.0251768272973862)
  566.   (f2cl-lib:fset (f2cl-lib:fref gama (16) ((1 26))) 0.023457075530607902)
  567.   (f2cl-lib:fset (f2cl-lib:fref gama (17) ((1 26))) 0.0219508390134907)
  568.   (f2cl-lib:fset (f2cl-lib:fref gama (18) ((1 26))) 0.020621082823564597)
  569.   (f2cl-lib:fset (f2cl-lib:fref gama (19) ((1 26))) 0.0194388240897881)
  570.   (f2cl-lib:fset (f2cl-lib:fref gama (20) ((1 26))) 0.018381063380068303)
  571.   (f2cl-lib:fset (f2cl-lib:fref gama (21) ((1 26))) 0.0174293213231963)
  572.   (f2cl-lib:fset (f2cl-lib:fref gama (22) ((1 26))) 0.0165685837786612)
  573.   (f2cl-lib:fset (f2cl-lib:fref gama (23) ((1 26))) 0.015786528598791802)
  574.   (f2cl-lib:fset (f2cl-lib:fref gama (24) ((1 26))) 0.015072950149409601)
  575.   (f2cl-lib:fset (f2cl-lib:fref gama (25) ((1 26))) 0.0144193250839955)
  576.   (f2cl-lib:fset (f2cl-lib:fref gama (26) ((1 26))) 0.0138184805735342)
  577.   (defun dasyjy (funjy x fnu flgjy in y wk iflw)
  578.     (declare (type (array double-float (*)) wk y)
  579.              (type f2cl-lib:integer4 iflw in)
  580.              (type double-float flgjy fnu x)
  581.              (type single-float funjy)
  582.              (type
  583.               (function
  584.                (double-float double-float double-float double-float
  585.                 double-float)
  586.                (values &rest t))
  587.               funjy))
  588.     (f2cl-lib:with-array-data (y-%data% y-%offset% y)
  589.       (declare (type f2cl-lib:integer4 y-%offset%)
  590.                (type (simple-array double-float (*)) y-%data%)
  591.                (ignorable y-%offset% y-%data%))
  592.       (f2cl-lib:with-array-data (wk-%data% wk-%offset% wk)
  593.         (declare (type f2cl-lib:integer4 wk-%offset%)
  594.                  (type (simple-array double-float (*)) wk-%data%)
  595.                  (ignorable wk-%offset% wk-%data%))
  596.         (prog ((cr (make-array 10 :element-type 'double-float))
  597.                (dr (make-array 10 :element-type 'double-float))
  598.                (kmax (make-array 5 :element-type 'f2cl-lib:integer4))
  599.                (upol (make-array 10 :element-type 'double-float)) (abw2 0.0)
  600.                (akm 0.0) (ap 0.0) (asum 0.0) (az 0.0) (bsum 0.0) (crz32 0.0)
  601.                (dfi 0.0) (elim 0.0) (fi 0.0) (fn 0.0) (fn2 0.0) (phi 0.0)
  602.                (rcz 0.0) (rden 0.0) (relb 0.0) (rfn2 0.0) (rtz 0.0) (rzden 0.0)
  603.                (sa 0.0) (sb 0.0) (suma 0.0) (sumb 0.0) (s1 0.0) (ta 0.0)
  604.                (tau 0.0) (tb 0.0) (tfn 0.0) (tol 0.0) (t2 0.0) (xx 0.0) (z 0.0)
  605.                (z32 0.0) (i 0) (j 0) (jn 0) (jr 0) (ju 0) (k 0) (kb 0)
  606.                (klast 0) (kp1 0) (ks 0) (ksp1 0) (kstemp 0) (l 0) (lr 0)
  607.                (lrp1 0) (iseta 0) (isetb 0))
  608.           (declare (type (array f2cl-lib:integer4 (5)) kmax)
  609.                    (type f2cl-lib:integer4 isetb iseta lrp1 lr l kstemp ksp1 ks
  610.                     kp1 klast kb k ju jr jn j i)
  611.                    (type (array double-float (10)) upol dr cr)
  612.                    (type double-float z32 z xx t2 tol tfn tb tau ta s1 sumb
  613.                     suma sb sa rzden rtz rfn2 relb rden rcz phi fn2 fn fi elim
  614.                     dfi crz32 bsum az asum ap akm abw2))
  615.           (setf ta (f2cl-lib:d1mach 3))
  616.           (setf tol (max ta 1.0000000000000002e-15))
  617.           (setf tb (f2cl-lib:d1mach 5))
  618.           (setf ju (f2cl-lib:i1mach 15))
  619.           (if (= flgjy 1.0) (go label6))
  620.           (setf jr (f2cl-lib:i1mach 14))
  621.           (setf elim (* -2.303 tb (f2cl-lib:int-add ju jr)))
  622.           (go label7)
  623.          label6
  624.           (setf elim (* -2.303 (+ (* tb ju) 3.0)))
  625.          label7
  626.           (setf fn fnu)
  627.           (setf iflw 0)
  628.           (f2cl-lib:fdo (jn 1 (f2cl-lib:int-add jn 1))
  629.                         ((> jn in) nil)
  630.             (tagbody
  631.               (setf xx (/ x fn))
  632.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (1) ((1 *)) wk-%offset%)
  633.                              (- 1.0 (* xx xx)))
  634.               (setf abw2
  635.                       (coerce
  636.                        (abs (f2cl-lib:fref wk-%data% (1) ((1 *)) wk-%offset%))
  637.                        'double-float))
  638.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (2) ((1 *)) wk-%offset%)
  639.                              (f2cl-lib:fsqrt abw2))
  640.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (7) ((1 *)) wk-%offset%)
  641.                              (expt fn con2))
  642.               (if (> abw2 0.2775) (go label80))
  643.               (setf sa 0.0)
  644.               (if (= abw2 0.0) (go label10))
  645.               (setf sa (/ tols (f2cl-lib:flog abw2)))
  646.              label10
  647.               (setf sb sa)
  648.               (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  649.                             ((> i 5) nil)
  650.                 (tagbody
  651.                   (setf akm (max sa 2.0))
  652.                   (f2cl-lib:fset (f2cl-lib:fref kmax (i) ((1 5)))
  653.                                  (f2cl-lib:int akm))
  654.                   (setf sa (+ sa sb))
  655.                  label20))
  656.               (setf kb (f2cl-lib:fref kmax (5) ((1 5))))
  657.               (setf klast (f2cl-lib:int-sub kb 1))
  658.               (setf sa (f2cl-lib:fref gama (kb) ((1 26))))
  659.               (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
  660.                             ((> k klast) nil)
  661.                 (tagbody
  662.                   (setf kb (f2cl-lib:int-sub kb 1))
  663.                   (setf sa
  664.                           (+
  665.                            (* sa
  666.                               (f2cl-lib:fref wk-%data%
  667.                                              (1)
  668.                                              ((1 *))
  669.                                              wk-%offset%))
  670.                            (f2cl-lib:fref gama (kb) ((1 26)))))
  671.                  label30))
  672.               (setf z (* (f2cl-lib:fref wk-%data% (1) ((1 *)) wk-%offset%) sa))
  673.               (setf az (coerce (abs z) 'double-float))
  674.               (setf rtz (f2cl-lib:fsqrt az))
  675.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (3) ((1 *)) wk-%offset%)
  676.                              (* con1 az rtz))
  677.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (4) ((1 *)) wk-%offset%)
  678.                              (*
  679.                               (f2cl-lib:fref wk-%data% (3) ((1 *)) wk-%offset%)
  680.                               fn))
  681.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (5) ((1 *)) wk-%offset%)
  682.                              (* rtz
  683.                                 (f2cl-lib:fref wk-%data%
  684.                                                (7)
  685.                                                ((1 *))
  686.                                                wk-%offset%)))
  687.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (6) ((1 *)) wk-%offset%)
  688.                              (*
  689.                               (-
  690.                                (f2cl-lib:fref wk-%data%
  691.                                               (5)
  692.                                               ((1 *))
  693.                                               wk-%offset%))
  694.                               (f2cl-lib:fref wk-%data%
  695.                                              (5)
  696.                                              ((1 *))
  697.                                              wk-%offset%)))
  698.               (if (<= z 0.0) (go label35))
  699.               (if (> (f2cl-lib:fref wk-%data% (4) ((1 *)) wk-%offset%) elim)
  700.                   (go label75))
  701.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (6) ((1 *)) wk-%offset%)
  702.                              (-
  703.                               (f2cl-lib:fref wk-%data%
  704.                                              (6)
  705.                                              ((1 *))
  706.                                              wk-%offset%)))
  707.              label35
  708.               (setf phi (f2cl-lib:fsqrt (f2cl-lib:fsqrt (+ sa sa sa sa))))
  709.               (setf kb (f2cl-lib:fref kmax (5) ((1 5))))
  710.               (setf klast (f2cl-lib:int-sub kb 1))
  711.               (setf sb (f2cl-lib:fref beta (kb 1) ((1 26) (1 5))))
  712.               (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
  713.                             ((> k klast) nil)
  714.                 (tagbody
  715.                   (setf kb (f2cl-lib:int-sub kb 1))
  716.                   (setf sb
  717.                           (+
  718.                            (* sb
  719.                               (f2cl-lib:fref wk-%data%
  720.                                              (1)
  721.                                              ((1 *))
  722.                                              wk-%offset%))
  723.                            (f2cl-lib:fref beta (kb 1) ((1 26) (1 5)))))
  724.                  label40))
  725.               (setf ksp1 1)
  726.               (setf fn2 (* fn fn))
  727.               (setf rfn2 (/ 1.0 fn2))
  728.               (setf rden 1.0)
  729.               (setf asum 1.0)
  730.               (setf relb (* tol (abs sb)))
  731.               (setf bsum sb)
  732.               (f2cl-lib:fdo (ks 1 (f2cl-lib:int-add ks 1))
  733.                             ((> ks 4) nil)
  734.                 (tagbody
  735.                   (setf ksp1 (f2cl-lib:int-add ksp1 1))
  736.                   (setf rden (* rden rfn2))
  737.                   (setf kstemp (f2cl-lib:int-sub 5 ks))
  738.                   (setf kb (f2cl-lib:fref kmax (kstemp) ((1 5))))
  739.                   (setf klast (f2cl-lib:int-sub kb 1))
  740.                   (setf sa (f2cl-lib:fref alfa (kb ks) ((1 26) (1 4))))
  741.                   (setf sb (f2cl-lib:fref beta (kb ksp1) ((1 26) (1 5))))
  742.                   (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
  743.                                 ((> k klast) nil)
  744.                     (tagbody
  745.                       (setf kb (f2cl-lib:int-sub kb 1))
  746.                       (setf sa
  747.                               (+
  748.                                (* sa
  749.                                   (f2cl-lib:fref wk-%data%
  750.                                                  (1)
  751.                                                  ((1 *))
  752.                                                  wk-%offset%))
  753.                                (f2cl-lib:fref alfa (kb ks) ((1 26) (1 4)))))
  754.                       (setf sb
  755.                               (+
  756.                                (* sb
  757.                                   (f2cl-lib:fref wk-%data%
  758.                                                  (1)
  759.                                                  ((1 *))
  760.                                                  wk-%offset%))
  761.                                (f2cl-lib:fref beta (kb ksp1) ((1 26) (1 5)))))
  762.                      label50))
  763.                   (setf ta (* sa rden))
  764.                   (setf tb (* sb rden))
  765.                   (setf asum (+ asum ta))
  766.                   (setf bsum (+ bsum tb))
  767.                   (if (and (<= (abs ta) tol) (<= (abs tb) relb)) (go label70))
  768.                  label60))
  769.              label70
  770.               (setf bsum
  771.                       (/ bsum
  772.                          (* fn
  773.                             (f2cl-lib:fref wk-%data%
  774.                                            (7)
  775.                                            ((1 *))
  776.                                            wk-%offset%))))
  777.               (go label160)
  778.              label75
  779.               (setf iflw 1)
  780.               (go end_label)
  781.              label80
  782.               (f2cl-lib:fset (f2cl-lib:fref upol (1) ((1 10))) 1.0)
  783.               (setf tau
  784.                       (/ 1.0
  785.                          (f2cl-lib:fref wk-%data% (2) ((1 *)) wk-%offset%)))
  786.               (setf t2
  787.                       (/ 1.0
  788.                          (f2cl-lib:fref wk-%data% (1) ((1 *)) wk-%offset%)))
  789.               (if (>= (f2cl-lib:fref wk-%data% (1) ((1 *)) wk-%offset%) 0.0)
  790.                   (go label90))
  791.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (3) ((1 *)) wk-%offset%)
  792.                              (coerce
  793.                               (abs
  794.                                (-
  795.                                 (f2cl-lib:fref wk-%data%
  796.                                                (2)
  797.                                                ((1 *))
  798.                                                wk-%offset%)
  799.                                 (atan
  800.                                  (f2cl-lib:fref wk-%data%
  801.                                                 (2)
  802.                                                 ((1 *))
  803.                                                 wk-%offset%))))
  804.                               'double-float))
  805.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (4) ((1 *)) wk-%offset%)
  806.                              (*
  807.                               (f2cl-lib:fref wk-%data% (3) ((1 *)) wk-%offset%)
  808.                               fn))
  809.               (setf rcz
  810.                       (/ (- con1)
  811.                          (f2cl-lib:fref wk-%data% (4) ((1 *)) wk-%offset%)))
  812.               (setf z32
  813.                       (* 1.5
  814.                          (f2cl-lib:fref wk-%data% (3) ((1 *)) wk-%offset%)))
  815.               (setf rtz (expt z32 con2))
  816.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (5) ((1 *)) wk-%offset%)
  817.                              (* rtz
  818.                                 (f2cl-lib:fref wk-%data%
  819.                                                (7)
  820.                                                ((1 *))
  821.                                                wk-%offset%)))
  822.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (6) ((1 *)) wk-%offset%)
  823.                              (*
  824.                               (-
  825.                                (f2cl-lib:fref wk-%data%
  826.                                               (5)
  827.                                               ((1 *))
  828.                                               wk-%offset%))
  829.                               (f2cl-lib:fref wk-%data%
  830.                                              (5)
  831.                                              ((1 *))
  832.                                              wk-%offset%)))
  833.               (go label100)
  834.              label90
  835.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (3) ((1 *)) wk-%offset%)
  836.                              (coerce
  837.                               (abs
  838.                                (-
  839.                                 (f2cl-lib:flog
  840.                                  (/
  841.                                   (+ 1.0
  842.                                      (f2cl-lib:fref wk-%data%
  843.                                                     (2)
  844.                                                     ((1 *))
  845.                                                     wk-%offset%))
  846.                                   xx))
  847.                                 (f2cl-lib:fref wk-%data%
  848.                                                (2)
  849.                                                ((1 *))
  850.                                                wk-%offset%)))
  851.                               'double-float))
  852.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (4) ((1 *)) wk-%offset%)
  853.                              (*
  854.                               (f2cl-lib:fref wk-%data% (3) ((1 *)) wk-%offset%)
  855.                               fn))
  856.               (setf rcz
  857.                       (/ con1
  858.                          (f2cl-lib:fref wk-%data% (4) ((1 *)) wk-%offset%)))
  859.               (if (> (f2cl-lib:fref wk-%data% (4) ((1 *)) wk-%offset%) elim)
  860.                   (go label75))
  861.               (setf z32
  862.                       (* 1.5
  863.                          (f2cl-lib:fref wk-%data% (3) ((1 *)) wk-%offset%)))
  864.               (setf rtz (expt z32 con2))
  865.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (7) ((1 *)) wk-%offset%)
  866.                              (expt fn con2))
  867.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (5) ((1 *)) wk-%offset%)
  868.                              (* rtz
  869.                                 (f2cl-lib:fref wk-%data%
  870.                                                (7)
  871.                                                ((1 *))
  872.                                                wk-%offset%)))
  873.               (f2cl-lib:fset (f2cl-lib:fref wk-%data% (6) ((1 *)) wk-%offset%)
  874.                              (*
  875.                               (f2cl-lib:fref wk-%data% (5) ((1 *)) wk-%offset%)
  876.                               (f2cl-lib:fref wk-%data%
  877.                                              (5)
  878.                                              ((1 *))
  879.                                              wk-%offset%)))
  880.              label100
  881.               (setf phi (f2cl-lib:fsqrt (* (+ rtz rtz) tau)))
  882.               (setf tb 1.0)
  883.               (setf asum 1.0)
  884.               (setf tfn (/ tau fn))
  885.               (setf rden (/ 1.0 fn))
  886.               (setf rfn2 (* rden rden))
  887.               (setf rden 1.0)
  888.               (f2cl-lib:fset (f2cl-lib:fref upol (2) ((1 10)))
  889.                              (*
  890.                               (+ (* (f2cl-lib:fref c (1) ((1 65))) t2)
  891.                                  (f2cl-lib:fref c (2) ((1 65))))
  892.                               tfn))
  893.               (setf crz32 (* con548 rcz))
  894.               (setf bsum (+ (f2cl-lib:fref upol (2) ((1 10))) crz32))
  895.               (setf relb (* tol (abs bsum)))
  896.               (setf ap tfn)
  897.               (setf ks 0)
  898.               (setf kp1 2)
  899.               (setf rzden rcz)
  900.               (setf l 2)
  901.               (setf iseta 0)
  902.               (setf isetb 0)
  903.               (f2cl-lib:fdo (lr 2 (f2cl-lib:int-add lr 2))
  904.                             ((> lr 8) nil)
  905.                 (tagbody
  906.                   (setf lrp1 (f2cl-lib:int-add lr 1))
  907.                   (f2cl-lib:fdo (k lr (f2cl-lib:int-add k 1))
  908.                                 ((> k lrp1) nil)
  909.                     (tagbody
  910.                       (setf ks (f2cl-lib:int-add ks 1))
  911.                       (setf kp1 (f2cl-lib:int-add kp1 1))
  912.                       (setf l (f2cl-lib:int-add l 1))
  913.                       (setf s1 (f2cl-lib:fref c (l) ((1 65))))
  914.                       (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
  915.                                     ((> j kp1) nil)
  916.                         (tagbody
  917.                           (setf l (f2cl-lib:int-add l 1))
  918.                           (setf s1
  919.                                   (+ (* s1 t2) (f2cl-lib:fref c (l) ((1 65)))))
  920.                          label110))
  921.                       (setf ap (* ap tfn))
  922.                       (f2cl-lib:fset (f2cl-lib:fref upol (kp1) ((1 10)))
  923.                                      (* ap s1))
  924.                       (f2cl-lib:fset (f2cl-lib:fref cr (ks) ((1 10)))
  925.                                      (* (f2cl-lib:fref br (ks) ((1 10)))
  926.                                         rzden))
  927.                       (setf rzden (* rzden rcz))
  928.                       (f2cl-lib:fset (f2cl-lib:fref dr (ks) ((1 10)))
  929.                                      (* (f2cl-lib:fref ar (ks) ((1 8))) rzden))
  930.                      label120))
  931.                   (setf suma (f2cl-lib:fref upol (lrp1) ((1 10))))
  932.                   (setf sumb
  933.                           (+
  934.                            (f2cl-lib:fref upol
  935.                                           ((f2cl-lib:int-add lr 2))
  936.                                           ((1 10)))
  937.                            (* (f2cl-lib:fref upol (lrp1) ((1 10))) crz32)))
  938.                   (setf ju lrp1)
  939.                   (f2cl-lib:fdo (jr 1 (f2cl-lib:int-add jr 1))
  940.                                 ((> jr lr) nil)
  941.                     (tagbody
  942.                       (setf ju (f2cl-lib:int-sub ju 1))
  943.                       (setf suma
  944.                               (+ suma
  945.                                  (* (f2cl-lib:fref cr (jr) ((1 10)))
  946.                                     (f2cl-lib:fref upol (ju) ((1 10))))))
  947.                       (setf sumb
  948.                               (+ sumb
  949.                                  (* (f2cl-lib:fref dr (jr) ((1 10)))
  950.                                     (f2cl-lib:fref upol (ju) ((1 10))))))
  951.                      label130))
  952.                   (setf rden (* rden rfn2))
  953.                   (setf tb (- tb))
  954.                   (if (> (f2cl-lib:fref wk-%data% (1) ((1 *)) wk-%offset%) 0.0)
  955.                       (setf tb (coerce (abs tb) 'double-float)))
  956.                   (if (< rden tol) (go label131))
  957.                   (setf asum (+ asum (* suma tb)))
  958.                   (setf bsum (+ bsum (* sumb tb)))
  959.                   (go label140)
  960.                  label131
  961.                   (if (= iseta 1) (go label132))
  962.                   (if (< (abs suma) tol) (setf iseta 1))
  963.                   (setf asum (+ asum (* suma tb)))
  964.                  label132
  965.                   (if (= isetb 1) (go label133))
  966.                   (if (< (abs sumb) relb) (setf isetb 1))
  967.                   (setf bsum (+ bsum (* sumb tb)))
  968.                  label133
  969.                   (if (and (= iseta 1) (= isetb 1)) (go label150))
  970.                  label140))
  971.              label150
  972.               (setf tb (f2cl-lib:fref wk-%data% (5) ((1 *)) wk-%offset%))
  973.               (if (> (f2cl-lib:fref wk-%data% (1) ((1 *)) wk-%offset%) 0.0)
  974.                   (setf tb (- tb)))
  975.               (setf bsum (/ bsum tb))
  976.              label160
  977.               (multiple-value-bind
  978.                   (var-0 var-1 var-2 var-3 var-4)
  979.                   (funcall funjy
  980.                            (f2cl-lib:fref wk-%data% (6) ((1 *)) wk-%offset%)
  981.                            (f2cl-lib:fref wk-%data% (5) ((1 *)) wk-%offset%)
  982.                            (f2cl-lib:fref wk-%data% (4) ((1 *)) wk-%offset%)
  983.                            fi
  984.                            dfi)
  985.                 (declare (ignore))
  986.                 (when var-0
  987.                   (f2cl-lib:fset
  988.                    (f2cl-lib:fref wk-%data% (6) ((1 *)) wk-%offset%)
  989.                    var-0))
  990.                 (when var-1
  991.                   (f2cl-lib:fset
  992.                    (f2cl-lib:fref wk-%data% (5) ((1 *)) wk-%offset%)
  993.                    var-1))
  994.                 (when var-2
  995.                   (f2cl-lib:fset
  996.                    (f2cl-lib:fref wk-%data% (4) ((1 *)) wk-%offset%)
  997.                    var-2))
  998.                 (when var-3 (setf fi var-3))
  999.                 (when var-4 (setf dfi var-4)))
  1000.               (setf ta (/ 1.0 tol))
  1001.               (setf tb (* (f2cl-lib:d1mach 1) ta 1000.0))
  1002.               (if (> (abs fi) tb) (go label165))
  1003.               (setf fi (* fi ta))
  1004.               (setf dfi (* dfi ta))
  1005.               (setf phi (* phi tol))
  1006.              label165
  1007.               (f2cl-lib:fset (f2cl-lib:fref y-%data% (jn) ((1 *)) y-%offset%)
  1008.                              (/ (* flgjy phi (+ (* fi asum) (* dfi bsum)))
  1009.                                 (f2cl-lib:fref wk-%data%
  1010.                                                (7)
  1011.                                                ((1 *))
  1012.                                                wk-%offset%)))
  1013.               (setf fn (- fn flgjy))
  1014.              label170))
  1015.           (go end_label)
  1016.          end_label
  1017.           (return (values nil nil nil nil nil nil nil iflw)))))))
  1018.  
  1019.